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H. Matsuzaki and K. Yamakoshi 
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To study origin and evolution of the interplanetary dust, it is very important to investigate the size 
distribution. Here the changes of the size distribution of meteoroid particles due to the ablative effects during 
atmospheric entry were investigated by numerical computer simulation. Using the results, the pre-atmosphcric 
size distribution of the interplanetary dust particles could be estimated from that of ablated spherules taken from 
deep-sea sediments. Wc arc now analyzing deep-sea spherules from some aspects and examining if we could get 
any information about the interplanetary dust. 

Fundamental equations of the numerical simulation are followed[ 1 ] ; 

( 1 ) dp/dt = Mdv/dt = -p atm vAv + Mg, (2) 1/2 p atm Av 3 = jis 2 e oT 1 - LdM/dt, 

(3) A = 1/4 Jts\ (4) M = 1/6 ;ips\ (5) -dM/dt = jis 2 p e (pm H /2jikT) 1 \ 

(P, M, v, s, p: momentum, mass, velocity, diameter, density of the meteoroid particle; p alm : atmospheric 
density; g: gravitational acceleration; f : cmissivily(=i.O); a: Stephen-Boltzmann constant; L, p e , p: latent heat 
of vaporization, vapor pressure, molecular weight of the meteoroid material; m H : mass of the hydrogen atom; k: 

Boltzmann coastant) 

(1) is the equation of motion, in which the momentum change of the meteoroid particle going through the 
atmosphere is equal to sum of the momenta of colliding air molecules and gravitational effect. Here, we can set 
dp=Mdv if evaporating mass is to dissipate in 4ji direction. (2) represents energy balance. Kinetic energy 
transported to the meteoroid particle by colliding air molecules is to be dissipated through radiation and ablative 
mass loss. We set (3) and (4) so that the maeteoroid particle is spherical. Mass loss rate is determined by (3). p c 
is a function of temperature of the meteoroid, T, such that; p e (T) = l(F Uc2 ' 1 (in MKS units). So, in practice, we 
solve (2) as the equation of T. The value of p at each altitude is taken from U. S. Standard Atmosphere of 
1976. In this study, two kinds of meteoroid particles, iron and stony, were taken. Physical properties of the 
materials used in this calculation are presented in Table 1 [2,3 J. The calculation starts at the meteoroid altitude 
190km at t=0, then equations (1) to (5) are solved at every 0.02 seconds. The calculation is terminated if the 
altitude is over 250km or below 50km, or the velocity becomes less than lkm/s. 

Parameters in this calculation are the initial size of the matcoroid particle, s 0 , atmospheric entry angle, 0 
(in the case of perpendicular input, 0=90°), and initial velocity, v Q . 7920 atmospheric entry events with various 
(s o ,0,v o ) value combinations ranged 10pm < s 0 < 1000pm, 0° < 0 < 90°, and 11.2km/s < v 0 < 20.2km/s were 
calculated. To these 7920 events, pre-atmosphcric (initial) size distribution; dn « s u A ds 0 , entry angle 
distribution; dn a - d(cos20)[ I ], and initial velocity distribution; dn « v 0 ‘ v \ivj4], were taken into account, 
where dn represents the number of particles (or probability) having initial size between s 0 and s +ds 0 , entry 
angle between 0 and 0+d0, and initial velocity between v y and v 0 +dv 0 , respectively. Index value A in the initial 
size distribution, we call "the initial index A n here after, was varied from 2 to 7. To each initial index A, the 
final size distribution was counted out in accordance with these three distributions. For example, the final size 
distribution was found to be shown as Fig. 1 in the case of iron material and the initial index A=4.5. (In this 
diagram, relative population before atmospheric entry is to be unity \v i t h parameter s u between 45 pm and 55 p- 
m, and parameter v 0 between 13.7km/s and I4.7km/s.) In Fig. 1 we can see that the final size distribution is 
also expressed nearly as; dn « s’Vs, and the final index B is found to be 5.9 in this case. Fig. 2 was obtained 
by calculations in cases of various initial indices A and counting out the correspoding final indices B. The 
relation between the initial index A and the final index B was found to be formally as; B = 1.20A + 0.40 (Iron). 

Similarly, about stony material particles, the relation was such that; B = 1.45A - 0.25 (Stony). 

Using these relations, we can estimate the pre-atmosphcric size distributions of iron and stony meteoroid 
particles inversely by measuring that of ablated iron and stony spherules taken from deep-sea sediments. Both 
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iron and stony spherules were picked up from the same fraction of sediments which had been taken as a dredged 
sample : GH79-1, Sl.1476 (9°49.96'N, 167°17.68’W, 5151m-dcpth). In this work, iron spherules in size larger 
than 80pm were counted, whereas many smaller ones exsisted. Stony spherules smaller than 100pm were rarely 
found, so the almose of all stony spherules could be thought to be picked up from the fraction. The results were 
showed in Fig. 3 in cumulative numbers with some other studies’ [5,6,7]. Iron spherules are found to have the 
gradient -1.9— 2.5 in Fig. 3, that is, the final index B = 2 .9-3. 5 (Iron). Using the relation oblined before, we 
get the initial index A = 2. 1-2.6 (Iron). Similarly, in the case of stony spherules, the final index is; B = 
4.8-5.0 (Stony), and then the initial index is; A = 3.5-3.6 (Stony). 

Dohnanyi (1972)[8] gave the size distribution of the interplanetary meteoroid particles as; dn s* 4 5 ds 
(s-lOOpm). More recently, we can get from the data compilation of Grun et al. (1985) [9]; dn s' 4 °~‘ 4 6 ds (s-100 
pm). These size distributions are consistent with that estimated from the stony spherules. However, the iron 
spherules themselves have smaller index of size distribution (2.9-3.5) so that estimated index of pre- 
atmospheric size distribution (2. 1-2.6) is much smaller than that of the interplanetary meteoroid particles 
(4.0-4 .6). As a consequence, it is suggested that the major components of the interplanetary meteoroid particles 
are stony particles, and that, in the interplanetary space, iron particles have different size distribution from stony 
ones in the size range around 100pm. 

(References) [1] Love, J. J. and Brownlee, D. E. (1991) Icarus, 89, 26-43; [2] Bronshten, A. E. (1983) 
"Physics of meteoric phenomena", Geophys. Astrophys. Monographs, D. Reidel; [3] Kornblum, J. J. (1969) 
JGR, 74(8 ), 1893-1907; [4] Fraundorf, P. (1980) Geophys. Res. Let., 10, 765-768; [5] Brownlow, A. E. et al. 
(1966) Geophys. J. R. Astro. Soc., 12, 1-12; [6] Murrell, M. T. et al. (1980) GCA, 44, 2067-74; [7] 
Yamakoshi, K. et al. (1981) JGR, 86(B4) , 3129-32; [8] Dohnanyi, J. S. (1972) Icarus, 17, 1-48; [9] Griin, E. 
et al. (1985) Icarus, 62, 244-272 


Table 1, Physical properties used in the calculation 
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Material Density p (g/cm 3 ) 
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Figure 2. Relation between the initial index A 
and the Final index B 



Figure 1. The final size dislribution of iron meteoroid 
particles in the case of the initial index A = 4.5 
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Figure 3. Cumulative size distributions 
of deep-sea spherules 






